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Abstract 

The production of vector bosons at present and future hadron colliders 
will provide a crucial test for QCD and Standard Model physics. In this 
paper we improve parton shower simulations of Drell-Yan processes by 
implementing matrix-element corrections to the initial-state radiation. We 
apply our work to the HERWIG Monte Carlo event generator and compare 
our phenomenological results with the ones obtained using the previous 
version of HERWIG, with resummed calculations which we match to the 
exact first-order perturbative result, and with recent Tevatron data. We 
also make some predictions for jet events at the LHC. 
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1 Introduction 



The production of vector bosons W^, and 7* [|I| in high energy hadronic colhsions 
is one of the most important processes that should be investigated in order to test the 
Standard Model of electroweak interactions and Quantum Chromodynamics. By mea- 
suring the rapidity or the mass distribution of the leptonic decay products one can also 
investigate the quark and antiquark distributions inside the colliding hadrons. For such 
processes, generated at the lowest order via a hard scattering like qq' ^ V, where q and 
q' have the same flavour for Z/'-f* production and different flavour for W production, 
higher order corrections due to multiple gluon emission in the initial-state radiation will 
play a crucial role. Many analyses have been devoted to the phenomenology of vector 
boson production, particularly to the differential distribution with respect to the trans- 
verse momentum qx of the produced vector boson. The approach of resummation of 
large logarithms of the ratio my / qr has been followed in many cases. This was origi- 
nally proposed by Dokshitzer, Dyakonov and Troyan (DDT) and then accomplished 
by Collins, Soper and Sterman (CSS) [Q, who performed the leading logarithmic resum- 
mation in the space of the impact parameter b, which is the Fourier space conjugate 
to qx- Ladinsky and Yuan implemented the CSS results numerically @]. Resummations 
of the initial-state multiple emission have been performed in |Q in the 6-space and more 
recently in |0-1O| in both the b- and the g^-space. In the resummation is also 

matched with the exact perturbative first-order result, which is important at high qjn- In 
the 6-space approach non-perturbative effects in the region of large values of b are taken 
into account via Gaussian functions in b, corresponding to a smearing of the transverse 
momentum distribution 0,^, which can also be directly implemented in g-p-space f^. 

Another possible approach to studying the phenomenology of vector bosons is to use 
Monte Carlo simulations of the initial-state parton shower. Standard parton showers 
n] , |l2| are performed in the leading-log approximation, therefore they are reliable only 



in the soft or coUinear region of the phase space, corresponding to low qx values for 
the produced vector boson. If we wish to study the high qj- region of the spectrum it 
is necessary to provide parton showers with matrix-element corrections. Refs. [|T3],|T1| 
implement matrix-element corrections to simulations of vector boson production in the 
PYTHIA Monte Carlo event generator and compare them with the results obtained at 
the Tevatron collider by the D0 collaboration. 

In this paper we reconsider this problem and apply matrix-element corrections to the 
initial-state radiation of the HERWIG parton shower, following the general prescription 
contained in as we already did for e^e~ annihilation [0, Deep Inelastic Scattering 
[pT| and top quark decays 



It is worth recalling that at present no Monte Carlo program including the full next- 
to-leading order (NLO) results exists, as it is not known how to set up a full NLO 
calculation in a probabilistic way. When providing parton showers with matrix-element 
corrections we still only get the leading-order normalization, because in the initial-state 
cascade we only include leading logs and not the full one-loop virtual contributions. 

In Section 2 we review the basis of the HERWIG parton shower algorithm for the 
initial-state radiation in hadronic collisions. In Sections 3 and 4 we discuss the hard and 
soft matrix-element corrections to vector boson production. In Section 5 we plot some 
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relevant phenomenological distributions at the centre-of-mass energy of the Tevatron 
and of the LHC using the new version of HERWIG. We compare our results with 
previous versions of HERWIG, resummed calculations and experimental data. Finally, 
in Section 6 we discuss our results and make some concluding comments. 

2 The parton shower algorithm 

The production of a vector boson V in hadronic collisions is given at lowest order by 
the elementary parton-level process qq V. In the following, we shall assume that the 
vector boson decays into a lepton pair (Drell-Yan interactions). The first-order tree- 
level corrections to such a process are given by the processes qq Vg and qg Vq 
{qg — > Vq), where the initial-state partons can come from either incoming hadron. 

A possible method to implement the initial-state parton shower in a probabilistic 
way is the Altarelli-Parisi approach, in which the initial energy scale Qo is increased up 
to the probed value Q and all the effect of the emitted partons is integrated out. 

On the contrary, standard Monte Carlo programs [|T^,0 explicitly keep track of the 
accompanying radiation, by implementing the so-called 'backward evolution' in which 
the hard scale is reduced away from the hard vertex, tracing the hard scattering partons 
back into the original incoming hadrons and explicitly generating the distribution of 
emitted partons. In the leading infrared approximation, the probability of the emission 
of an additional parton from a parton i is given by the general result for the radiation 
of a soft /collinear parton: 



The HERWIG parton shower is ordered according to the variable qf = E'^^i, where 
E is the energy of the parton that split and = where Pi is the four-momentum 
of the emitted parton; ph is a lightlike vector with momentum component parallel to 
the incoming hadron; Eh and Ei are the energy components of ph and pf, and Zi is 
the energy fraction of the outgoing space-like parton (which goes on to participate 
in the hard process) with respect to the incoming one (i.e. Zi = 1 — E.i/E). In the 
approximation of massless partons, we have = 1 — cos^, where 9 is the emission angle 
from the incoming hadron direction. When all emission is soft, the energy of the emitted 
partons is negligible [Ei -C E) , therefore ordering according to corresponds to angular 
ordering; when the emission is hard, the energy of the radiated parton is similar to that 
of the splitting parton, so ordering is equivalent to transverse momentum ordering. 
In (|l|) fa{xi, qf) is the parton distribution function for the partons of type a in the 
initial-state hadron, Xj being the parton energy fraction. At each step of the backward 
evolution a parton of type a, a quark for example, can evolve back to any other type 
of parton b, in this case either a quark of the same flavour or a gluon, having a higher 
value of Xi. 

The quantity /S.s,a{,q^,qV) the Sudakov form factor, resulting from the leading- 
logarithmic resummation and representing the probability that no resolvable radiation 
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is emitted from a parton of type a whose upper limit on emission is q^, with being, in 
the case of HERWIG, a cutoff on transverse momentum. This cutoff imphes a minimum 
value of the evolution scale q]; that can be reached, q^ > 4g^, but in practice this is 
smaller than the smallest scale at which most standard parton distribution function 
sets are reliable, so an additional cutoff on qf has to be applied. The ratio of form 
factors appearing in Eq. (|I|) represents the probability that the emission considered is 
the first, i.e. the one with the highest value of qf. In terms of Feynman diagrams, the 
Sudakov form factor sums up all-order virtual and unresolved contributions. Pabi^i) is 
the Altarelli-Parisi splitting function for a parton of type b to evolve to one of type a 
with momentum fraction Zi. as is the strong coupling, evaluated at a scale of order the 
transverse momentum of the emitted parton, which sums large higher order corrections. 
This, together with the angular ordering condition, makes HERWIG accurate to next- 
to-leading order at large x ||21|| . 

The definition of the variables qf and Zi is not Lorentz-invariant, but it is frame- 
dependent. Colour coherence implies that for any pair of colour-connected partons i and 
j the maximum values of the q variables are related by gimaxQ'jmax = Pi'Pj- Therefore one 
is free to choose the frame in which to define the initiating values g^max and g^max, with the 
only prescription being that their product must equal Pi ■ pj. The subsequent emissions 
are then ordered in qf. For vector boson production, as in most cases, symmetric limits 
are fixed by HERWIG, i.e. qf^^y, = qj^g^^ = Pi ■ Pj and the energy of the parton which 
initiates the cascade is set to E = gmax = y/Pi ■ Pj ■ Ordering according to qf therefore 
dictates < zf. 

After we generate the initial-state shower, the original partons are not on their 
mass-shell anymore, so their energy and momentum cannot be conserved. Energy- 
momentum conservation is then achieved by applying a separate boost to each jet along 
its own direction. As a result of this, the jet momenta are no longer equal to the parton 
ones, but energy and momentum are globally conserved and the vector boson acquires a 
transverse momentum from the recoil against the emitted partons. Since the mass shift 
becomes negligible in the soft and collinear limits, the precise details of this kinematic 
reshuffling are not fixed a priori, but are free choices of the model. 

Once the backward evolution has terminated, a model to reconstruct the original 
hadron is required. In HERWIG, if the backward evolution has not resulted in a va- 
lence quark, additional non-perturbative parton emission is generated to evolve back to 
a valence quark. Such a valence quark has a Gaussian distribution with respect to the 
non-perturbative intrinsic transverse momentum in the hadron, with a width that is an 
adjustable parameter of the model. In the following, when discussing the phenomeno- 
logical implications of our work, we shall consider both HERWIG's default value of zero, 
and an increased value of 1 GeV, bracketing the reasonable range of non-perturbative 
effects. 

The algorithm so far discussed is reliable only in the soft or collinear limits and, 
since it only describes radiation for < zf, there are regions of the phase space that 
are completely empty ('dead zones'). The radiation in such regions, according to the 
full matrix element, should be suppressed, but not completely absent as happens in 
HERWIG. We therefore need to improve the HERWIG model by implementing matrix- 
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element corrections. As usual jlS HT^] , this method works in two steps: we populate the 
missing phase space region by generating radiation according to a distribution obtained 
from the first-order matrix-element calculation ('hard corrections'); we correct the algo- 
rithm in the already-populated region using the matrix-element probability distribution 
whenever an emission is capable of being the 'hardest so far' ('soft corrections'). 



3 Hard Corrections 

In order to implement the hard and soft matrix-element corrections to simulations of 
the initial-state radiation in Drell-Yan processes, we firstly have to relate the HERWIG 
variables ^ and z to the kinematic ones we use in the matrix-element calculation. For 
the process q{pi)q{p2) — ^ giPsWi^) we parametrize the phase space according to the 
Mandelstam variables s = {pi + ^2)^, t = {pi — PsY and u = {p2 — PsY- Throughout 
this paper, we neglect the parton masses, so we have s + t + u = my. 
The phase space limits, in terms of the variables s and t, are: 

rriy < s < s, (2) 
my - s < i < 0, (3) 

where s is the squared energy in the centre-of-mass frame. Note that the point s = my 
corresponds to the soft singularity, and the lines t = and t = my — s to collinear 
emission. 

In order to relate s and t to ^ and z we use the property that the mass m and the 
transverse momentum pt of the q-g (q-g) jets are conserved in the showering frame. In 
doing this, we observe that, in the approximation of massless partons, the energy of 
the annihilating qq pair which produces the vector boson V is equal to E' = ^Pq ■ Pg = 



m^ = —^ml,, (4) 



In terms of the showering variables, we obtain: 

1 

z^ 

V\ = ^^^e(2-0<, (5) 
and in terms of the matrix-element variables: 

= i, (6) 

= T (7) 

Combining them we get the following equations: 



m 



ml\ ^ 



i ^ i ) si 



""v _ '"•v"" _|_ rn'"^ /_L _ ^ J 

m?_ r, rn n — * V / 



4 




2^ 

rriv 



Figure 1: The phase space for additional emission in Drell-Yan production with s/my = 
6.25, showing the kinematic hmits (sohd) and HERWIG's parton shower hmits (dashed). 

The region where HERWIG does not allow gluon radiation can be derived by solving 
the equation ^ > z^: 

Smin < S < S (10) 
^min ^ t irnax) (-^-^) 

where tmax can be obtained by solving the equation 

P + 3mli + 2m\r (l-'^\= 0, (12) 



or: 



tmax = -^ ( 3 - + 8m2./s ) . (13) 



It is straightforward to write tmin as 



while Smin can be determined by the condition tmin('S) < ^max( 

m 
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^ (7 - . (15) 



In Fig. m we plot the total phase space and HERWIG's limits for a vector boson 
mass of my = 80 GeV and a centre-of-mass energy of -y/i = 200 GeV, in terms of 



the normalized variables s/rriy and t/rriy. We can see that, as in cases |jT6| and JT? 



and differently from |18[, the soft and the coUinear singularities are well inside the 



HERWIG phase space: we also have an overlapping region, corresponding to a kinematic 
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configuration in s and i where radiation can come from eitlier parton. Note tliat tlie 
only dependence on tlie external pliysical parameters is tlie position of tlie edge at large 
s, i.e. Smax = s, wfiile tlie values of Smin/my and the limits in t/rriy are independent 
of the centre-of-mass energy and of other kinematic conditions like the vector boson 
rapidity. 

Once we have the total and HERWIG phase space limits, in order to implement 
matrix-element corrections, we have to apply the exact differential cross section in the 
dead zone. In a general prescription is given to allow first-order corrections to 



quark scattering and annihilation processes once a generator of the lowest order process 
is available. For the Drell-Yan case, assuming that the virtuality and the rapidity of the 
produced vector boson are fixed by the Born process qq —>■ V, the first-order differential 
cross section da is proportional to the parton-level lowest order (Jq according to the 
relation: 

d^a = , /.A(^i)/.72fe)C.a, |cg J _ ^ _ _ 
fq/i{r]i)fq/2i^V 27r sHu L J 

In the above equation fq/i{xi) and fqiiixi) are the parton distribution functions of the 
scattering partons inside the incoming hadrons 1 and 2 for energy fractions Xi and X2 
in the process qq — > Vg^ while jqi\{j]\) and fqiiij]'!) refer to the Born process and cancel 
off the factors that are already in a^. The assumption that the rapidity and mass of the 
vector boson are the same as in the process qq allows us to recover the Born result 
in the limit of an extremely soft gluon radiation. 

As stated here, Eq. ([T6|) is a trivial rewriting of the first-order differential cross 
section, but the main point of p2| is that if the azimuth of the emitted gluon is generated 
in the right way, Eq. (|I6|) correctly describes the full process including the vector boson 
decay. Thus, to implement our matrix-element corrections, we do not need to know 
anything about the final state of the vector boson - its properties are correctly inherited 
from the Born process. 

In a similar way we deal with the Compton process q{pi)g{p3) — > '?(P2)^(q')- We 
define the Mandelstam variables s = {pi + Ps)^, t = {ps — P2Y and u = {pi — ^2)^, and 
we find the same expressions for the variables ^ and z and for the phase space limits. 
We obtain for the differential cross section [p2| : 

72 fq/l{Xl)fg/2{X2)TRas ds di r 2 A2 , / 2 -^^2] 

fq/l{Vl)fq/2{V2) 2n sH ^ J 

Extending this formula to processes where we have an incoming antiquark or where 
the gluon belongs to hadron 1 is straightforward. We then generate events according to 
the above distributions in the dead zone using standard techniques. 

When applying the hard corrections, in principle one should also implement the form 
factor, but, since we are quite far from the soft and collinear singularities, it is actually 
not important and we shall neglect it in the following. This is justified by the fact that 
the total fraction of events that receive an emission from the hard correction is small. 
For example for W production it is 3.9% at the Tevatron and 9.2% at the LHC. Also, 



6 



the fact that such fractions are quite small allows us to neglect multiple emissions in the 
dead zone and makes the use of the exact first-order result reliable. 

Among these events, the fraction of qq' — > Wg processes is 53.5% at the Tevatron 
and 24.5% at the LHC The reason for the differences between the two machines can be 
understood in terms of the parton distribution functions. The gluon density inside the 
protons is higher when the colliding energy is increased because x is decreased; moreover 
at the LHC we have pp interactions instead of pp, therefore a qq' annihilation requires 
an antiquark q' to be taken from the 'sea'. 

The equivalent numbers for Z production are essentially identical. 



4 Soft Corrections 



According to [15|, we should also correct the emission in the region that is already 
populated by HERWIG using the exact first-order calculation for every emission that is 
the hardest so far. This can be performed by multiplying the parton shower distribution 
by a factor that is equal to the ratio of HERWIG's differential distribution to the matrix- 
element one. The only non-trivial part of this is in calculating the Jacobian factor 
J{s,i; z,^) of the transformation (z,^) {s,t). HERWIG's cross section is then given 
by 

(fa (fa - ^, 

where (fa/dzd^ is given by the elementary emission probability given in Eq. (p. The 
Jacobian factor J can be simply calculated from the relations given earlier: 

^ ti^TfijY ^) Z'^ 

J{s, t; z, ^) = — T77^ \2~/ I /-/I vT • V-'-^) 

s^ my^[l — z)^ [z + t,[l — z)) 

At this point we are able to make some comparisons with the approach that is 
followed in |jl3| where matrix-element corrections are added to the PYTHIA simulation 



of vector boson production. The parton shower probability distribution is applied over 
the whole phase space (in its older versions PYTHIA had a cutoff on the virtuality k"^ of 
the hard scattering parton that was constrained to be k"^ < in order to avoid double 
counting) and the exact matrix-element correction is applied only to the branching that 



is closest to the hard vertex. Unlike |jT3[, we have complementary phase space regions 
where we apply either the parton shower distribution (|l|) or the exact matrix-element 
ones ( P^ , p7| ), while in the parton shower region (^ < z^) we use the exact amplitude to 
generate the hardest emission so far instead of just the first emission. Correcting only 
the first emission can lead to problems due to the implementation of the Sudakov form 
factor whenever a subsequent harder emission occurs, as we would get the unphysical 
result that the probability of the hard emission would depend on the infrared cutoff that 
appears in the expression of the form factor. See |TH] for more details on this point. 
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5 Results 



Having implemented the hard and soft matrix-element corrections to the initial-state 
radiation in Drell-Yan interactions, we wish to investigate the impact they have on 
relevant phenomenological observables that can be measured at the Tevatron and in 
future at the LHC In the following, we shall mostly concentrate on W production 
although 7* and Z events are treated in exactly the same way. 



5.1 Vector boson transverse momentum 

A particularly significant phenomenological quantity is the transverse momentum of the 
W to the beam axis, which has been object of many theoretical and experimental anal- 
yses. In the soft/coUinear limit, the transverse momentum of the W is constrained to 
be qt < raw, since in the hard process qq' W the W is produced with no trans- 
verse momentum and it can acquire some qx only as a result of the initial-state parton 
showering. When the emission is generated according to the exact matrix element of 
processes like qq' —>■ Wg, qg q'W or qg — > q'W , the W produced in the hard process 
is allowed to have a non-zero qx and events with q^ > my/ are expected. 

In Fig. we compare the differential cross sections with respect to the W qr for 
pp collisions at the Tevatron energy[|, y/s = 1.8 TeV, obtained using HERWIG 5.9, the 
latest public version, with 6.1, the new version in progress where we include for the 
first time matrix-element corrections to the initial-state parton shower in vector boson 
production. We set an intrinsic transverse momentum equal to zero and use the MRS 



(Rl) parton distribution functions |2^. We can see from the plots that the impact of 
matrix-element corrections is negligible at low values of qx but it is quite relevant at 
high qx, where we have many more events with respect to the 5.9 version. Above some 
value of qT HERWIG 5.9 does not generate events anymore, while the 6.1 version still 
gives a non-zero differential cross section thanks to the events generated via the exact 



hard matrix element. As in e^e~ annihilation [jl6| and DIS[|T7l, it is actually the hard 
corrections that have a marked impact on our distributions, while the effect of the soft 
ones is quite negligible. 

It is also interesting to plot the qx spectrum obtained running the W + jets' process 
of HERWIG forcing the produced W to decay leptonically. This generates the hard 
process qq' Wg (or the equivalent ones with an initial-state gluon) for all events. As 
this matrix element diverges when the transverse momentum of the W approaches zero, 
HERWIG applies a user-defined cut on the generated in the hard process, which we 
set to 10 GeV. In our plot, this does not appear as a sharp cutoff since the W gets 
some recoil momentum due to the initial-state parton shower which can increase or 
decrease its transverse momentum. If, on the contrary, we had plotted the gr of the W 
generated in the hard 2 — > 2 process, we would have got a sharp peak at = 10 GeV. 
The agreement we find between the simulations of Drell-Yan processes provided with 
matrix-element corrections and the W + jet' events for large qr reassures us that the 



*The next Tevatron run will be at the slightly higher energy of 2 TeV. For the sake of comparison 
with existing data we use 1.8 TeV, but the results would not be qualitatively different for 2 TeV. 
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Figure 2: The W qx distribution at the Tevatron, according to HERWIG without (dot- 
ted) and with (sohd) matrix-element corrections. Also shown (dashed) is the W+ jet' 
process with a cutoff of 10 GeV. 



implementation of the hard corrections is reliable. 

We also see that the HERWIG distributions for Drell-Yan processes show a sharp 
peak in the first bin, which includes the value qt = 0: it corresponds to a fraction of 
events with no initial-state radiation and so W bosons produced with zero transverse 
momentum. With any fixed infrared cutoff value, one expects a non-zero (though ex- 
ponentially suppressed) fraction of events to give no resolvable radiation. These would 
normally be smeared out by non-perturbative effects like the intrinsic transverse mo- 
mentum, as we shall see in later plots, but since we set the width of its distribution to 
zero by default, all such events appear in the lowest qt bin. This is actually a technical 
deficiency of the Monte Carlo simulation and not a detectable physical effect. 

The D0 collaboration recently published data on the transverse momentum distri- 
bution of W bosons at the Tevatronp^. In Fig. ^ we compare the HERWIG predictions 
with it. In order to contribute to the investigation of possible effects of non-perturbative 
physics, we also run HERWIG setting an intrinsic partonic transverse momentum equal 
to 1 GeV, which we consider to be the maximum reasonable value. 

We see that the data has a significantly broader distribution at small qx than HER- 
WIG without intrinsic transverse momentum and that increasing the r.m.s. pt to 1 GeV 
is nowhere near enough to account for this. Furthermore the description of the data 
in the intermediate qt range, 30-70 GeV, is also rather poor. The intrinsic pt does 
not affect the predictions significantly for qt values above about 10 GeV, so there is no 
obvious way to improve the fit for intermediate values. 

However, these effects are actually because the D0 data are uncorrected for detector 
effects. We have run HERWIG through D0's fast simulation program, CMSp5[, and 
show results in Fig. ^. We see that HERWIG now describes the data rather well. The 
detector smearing is so strong at low qx that the additional smearing produced by the 
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Figure 3: The W qt distribution data from D0, in comparison with HERWIG with 
matrix-element corrections but without detector corrections. The sohd hne has zero 
intrinsic transverse momentum while the dashed one has an r.m.s. pt of 1 GeV. 
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Figure 4: As Fig. H but with the HERWIG results corrected for detector effects. 



intrinsic transverse momentum becomes irrelevant. 

At this point it is worthwhile commenting on the results shown in Refs. fl^jlj] • Both 
compare generator-level results with the D0 data. Ref. [Q's actually look rather similar 
to our generator-level results, so it is likely that after applying detector corrections they 
will describe the data as well as HERWIG. Ref. found good agreement with the D0 
data, but only after increasing the intrinsic transverse momentum to 4 GeV. It seems 
likely that this accounts for the smearing at low qt and would not be necessary after 



including detector smearing. However, in the intermediate qt range the results of |13 
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Figure 5: The W qt distribution at the LHC, according to HERWIG without (dotted) 
and with (sohd and dashed) matrix-element corrections, with zero intrinsic pt (sohd) 
and an r.m.s. pt of 1 GeV (dashed). 



are significantly lower than HERWIG. Since we find a detector correction of around 
a factor of two in this region, it would be very interesting to see the results of [O at 
detector level to see whether they are still able to fit the data. 

In Fig. 1^ we show the qx distributions for pp collisions at the energy of the LHC, 
^/s = 14 TeV, and find that the impact of the corrections is even bigger once the energy 
is increased. Unlike in the Tevatron transverse momentum distributions, we do not 
have the previously-mentioned sharp peak at qt = 0: this is because at the LHC we 
have pp interactions and the protons do not have valence antiquarks, while in order to 
produce a ly we do need a qq' hard scattering. As a result, the backward evolution 
has to produce at least one splitting, which always gives the W itself some transverse 
momentum. From the window in the top-right corner, we also see that at very low q^ 
the uncorrected version, 5.9, has a few percent more events than 6.1, particularly in 
the case of the LHC. As we said in the introduction, although we have matched to the 
tree-level NLO matrix elements, we still get the LO normalization, therefore the total 
cross sections obtained from versions 5.9 and 6.1 are the same. Since at the energy of 
the LHC we are generating a higher fraction of events at large g-r via the exact matrix- 
element distribution, it is reasonable that this enhancement is partially compensated by 
a slight suppression in the low qx region. 

In the region of low qT, it is worthwhile comparing the HERWIG distributions with 
some resummed calculations that are available in the literature. All these calcula- 
tions are based on the approach suggested in where the differential cross section 
with respect to the vector boson qT is expressed as the resummation of logarithms 
/ = log(my/g|.) to all orders in as- Two conflicting nomenclatures are used in the 
literature to denote which logarithms are summed: in the differential cross section, at 
each order in ag the largest term is ~ l/gf-Q;" which are sometimes known as the 
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Figure 6: The W qt distribution at the Tevatron, according to HERWIG with matrix- 
element corrections, with zero intrinsic pt (sohd histogram) and an r.m.s. pt of 1 GeV 
(dashed histogram), compared with the resummed results of [H in g^-space (solid) and 
in 6-space (dotted) and of in g^-space (dashed). 



leading logarithms, ~ l/gf. a" /2n-2 i^gj^g known as the next-to- leading logarithms, and 
so on. According to this classification, the results in ^ and |T0| are NNLL and NNNLL 
respectively. 

However, these logarithms 'exponentiate', allowing the differential cross section to be 
written in terms of the exponential (the 'form factor') of a series in as whose largest term 
is ~ a"/"'"'"^, which are also sometimes known as the leading logarithms, ~ a^/" being 
known as the next-to-leading logarithms, and so on. In |^ all NLL terms according to 
this nomenclature are summed in the form factor, which is evaluated either in the impact 
parameter 6-space or in the gy-space. The non-perturbative contribution is taken into 
account in the fe-space formalism following the general ideas in and setting Gaussian 
functions in the impact parameter b to quantify these effects. 

In Fig. 1^ we compare the HERWIG 6.1 differential cross section^ with the ones 
obtained from these approaches all normahzed to the corresponding total cross 

section. HERWIG clearly lies well within the range of the resummed approaches, except 
at very small qt where they become unreliable. To be more precise, the agreement is 
better between HERWIG and the two resummations in the g^-space, but even the b- 
space result is not too far from the Monte Carlo distribution. 

If we now wish to compare the HERWIG simulation with matrix-element corrections 

^As the resummed calculations deal with a fixed value of mw, Fig. ^is obtained running HERWIG 
with a vanishingly small W width, so mw — 80.4 GeV, its default value. W width effects are never- 
theless fully included in HERWIG and in the other plots we show. At low qr this assumption does not 
change the results dramatically, at high qt the effect of the W width is important as it allows values 
of qt larger that the default W mass even in the parton shower approximation, otherwise they could 
come only via the exact matrix-element generated events. 
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Figure 7: As Fig. ||, but with the resummed results matched with the exact 0{as) 
result. 

to the resummed calculations even for larger values of qt we need to match the latter 
with the exact 0{as) results to make them reliable there. This has already been done 
in the literature within the approach of |0J^, while in the analysis is limited to the 
low qt regime. 

Many prescriptions exist concerning how to perform such a matching. Ours is to sim- 
ply add the exact matrix-element cross sections for the parton level processes qq' —>■ Wg 
and q{q)g Wq'{q'), already calculated in (|TB|) and (|T^), to the resummed expressions 
and, in order to avoid double counting, subtract off the terms they have in common. It 
is straightforward to show that these are simply those terms in the exact 0{as) result 
that do not vanish as — 0. This prescription works fine if the resulting distribution 
is continuous at the point qx = mw, which means that the resummation and the low q^ 
0{ag) result exactly compensate each other and only the exact 'hard' matrix-element 
contribution survives. 

As discussed in it is not trivial to implement such a matching: the authors in 
fact do not succeed in obtaining a continuous distribution at the crucial matching point 
qx = rriw, but rather a step of size ~ a| was found. This comes about because the 
derivative of the Sudakov exponent is not required to go smoothly to zero at that point. 

We independently implement the matching for all the resummed calculations with 
which we wish to compare the HERWIG results and we find that the matching works well 
only for the g^-space resummation performed in [P]. For the 6-space and the approaches 
in P we do indeed find a step at qx = mw In fact even for the g^-space method of P 
we find a 'kink' at qx = i.e. although the curve is continuous, its derivative changes 
discontinuously there, albeit by an amount that is too small to notice on the figure. 

In Fig. 1^ we compare the HERWIG 6.1 distributions with 10-0 after the matching 
over the whole qx spectrum. We see that for the resummed distribution that is well 
matched and continuous for q^ = my/ the agreement with HERWIG is pretty good 
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Figure 8: The preliminary Z qx distribution data from CDF, in comparison with HER- 
WIG without (dotted) and with (sohd and dashed) matrix-element corrections. The 
solid and dotted lines have zero intrinsic transverse momentum while the dashed one 
has an r.m.s. pt of 1 GeV. 



everywhere. We do have slight discrepancies for medium values of g^; but they are well 
within the range that could be expected from the differences between the approaches 
followed by the Monte Carlo program and the calculation which keeps all the next-to- 
leading logarithms in the form factor. 

While the plots shown so far refer to W production, it is also worth comparing 
with some recent preliminary CDF data on Z production In Fig. ^, we have the 
CDF distribution with respect to the transverse momentum of the 7*/2' boson pro- 
duced at the Tevatron and decaying into an e^e~ pair with invariant mass in the range 
66 GeV < mz < 116 GeV. We compare the data with HERWIG before and after matrix- 
element corrections; we also normalize the HERWIG distribution to the experimental 
value of the cross section, 245.3 pb. The result is that we obtain good agreement with 
the experimental data only thanks to the application of the hard and soft corrections, 
otherwise the predictions would have been badly wrong for values qx > 50 GeV. There 
is perhaps some evidence that HERWIG does not produce enough smearing at low g^, 
even with an intrinsic pt of 1 GeV, with HERWIG peaking at about 2 GeV and the data 
peaking at about 3 GeV, but the overall fit is nevertheless acceptable. 

We have however found that better agreement can be obtained with an intrinsic pt 
of 2 GeV. 



5.2 Jet distributions 

We now look at the impact the matrix-element corrections have on the jet activity at 
the Tevatron and at the LHC. An interesting object to analyse is the hardest jet in 
transverse energy (the so-called 'first jet'). In Fig. ^ we plot the differential spectrum 
for the transverse energy of the first jet for ^/s = 1.8 TeV and y/s = 14 TeV, using 
HERWIG 5.9 and 6.1 and running the inclusive version of the fcy algorithm [0,^] for 
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Figure 9: The distribution in transverse momentum of the hardest jet in W production 
at the Tevatron (a) and the LHC (b) according to HERWIG without (dotted) and with 
(sohd) matrix-element corrections. 

a radius R = 0.5 at the Tevatron and R = 1 for the LHC |. The result is that the 
improvement introduced does have a significant effect since the number of events in 
which the first jet has high Et is significantly increased. The 5.9 and 6.1 distributions 
are similar for small values of Et, but for increasing Et the effect of the corrections 
introduced in HERWIG gets more and more relevant and at very high Et only events 
generated via the exact hard amplitude survive. The impact is really enormous in the 
case of the LHC as can be seen from Fig. ^d. 

In Fig. 1^ we plot the inclusive number of jets rijets that pass a transverse energy cut 
Et > 10 GeV. We see that implementing the matrix-element corrections significantly 
shifts the distribution towards larger rijcts- If we look at events with three or four jets 
having Et > 10 GeV, we see that their number is increased considerably both at the 
Tevatron and at the LHC. We have roughly an enhancement of a factor of 2 for three-jet 
events at both the energies, while for events with four high transverse energy jets, at 
the LHC we still get an enhancement of 2, while at the Tevatron the difference is almost 
a factor of 4. 



5.3 Rapidity distributions 

Fig. |ll] shows the distribution of the rapidity y of the dilepton pair at ^/s = 1.8 TeV 
and y/s = 14 TeV. As we sum over and W~ they are symmetric in ±y. We see that 
the matrix-element corrections do not significantly affect the rapidity distributions. 

Fig. shows the comparison between HERWIG and CDF for the rapidity of 
the produced e^e~ pair; the agreement is again good and the contribution of the matrix- 
element corrections is insignificant. 



"'"The correspondence bet'ween the radius R in the /ct algorithm and i?conc in an iterative cone 
algorithm is i?conc ~ 0.75 x R [ p8[ . The Tevatron experimentalists run an iterative cone algorithm 
■with radius i?cone = 0.4, so we choose R — 0.5 for the radius parameter "when "we consider jet events at 
^/s — 1.8 TeV. For the LHC "we stick to the recommended value of i? = 1. 
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Figure 10: The number of jets with Et > 10 GeV in W events at the Tevatron (a) and 
the LHC (b) according to HERWIG without (dotted) and with (sohd) matrix-element 
corrections. 
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Figure 11: The rapidity distribution of W bosons at the Tevatron (a) and the LHC (b) 
according to HERWIG without (dotted) and with (sohd) matrix-element corrections. 



6 Conclusions 

We have analysed Drell-Yan processes in hadron collisions in the Monte Carlo parton 
shower approach. This is accurate in the soft /coUinear approximation, but leaves empty 
regions in phase space. We implemented matrix-element corrections generating radiation 
according to the first-order amplitude in the dead zone and for every hardest so far emis- 
sion in the already populated region in the HERWIG parton shower. We compared our 
results with the previous version HERWIG 5.9, with experimental Tevatron data from 
the D0 and CDF collaborations and existing rcsummcd calculations of the spectrum 
of the transverse momentum of the vector boson. We found that the implemented 
corrections have a marked impact on the phenomenological distributions for high values 
of qt and the new version of HERWIG fits the D0 data for W production well over the 
whole qr spectrum after we correct the HERWIG results to detector level. At large qr 
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Figure 12: The Z rapidity distribution data from CDF, in comparison with HERWIG 
without (dotted) and with (sohd) matrix-element corrections. 

it is crucial to provide the Monte Carlo algorithm with matrix-element corrections in 
order to succeed in obtaining such an agreement. 

We also compared the HERWIG results after matrix-element corrections to some 
existing calculations based on a resummation of the initial-state radiation in the qT- 
space and in the 6-space. We found that in the range of low qT-, where actually the 
effect of matrix-element corrections is not so relevant and the resummed calculations 
are quite reliable, the parton shower distribution is in reasonable agreement with all of 
them, with discrepancies due to the methods followed by these different approaches. We 
also matched the resummed results to the exact Oipis) result, in such a way to allow 
them to be trustworthy at all qT values and found that the matching works well for a 
resummation performed in the q'r-space keeping all the next-to-leading logarithms in the 
Sudakov exponent. In this case, we also obtained good agreement with the HERWIG 6.1 

distribution. The other approaches considered showed a discontinuity at the point 
qT = mw once we match them to the exact first-order perturbative result. 

We also studied W + jet events at the Tevatron and at the LHC and found a 
significant effect of the new improvement of HERWIG as a larger number of jets of large 
transverse energy passes the typical experimental cuts. 

We compared the new version of HERWIG with the experimental data of the CDF 
collaboration on the transverse momentum and rapidity of Z bosons. We found good 
agreement after implementing the corrections. As a result, we feel confident that the 
simulation of vector boson production is now reliable. Using the new version of HERWIG 
6.1 to fit the experimental data will therefore provide us with better tests of the Standard 
Model and of QCD for the following Run II at the Tevatron and, ultimately, at the LHC. 

For the sake of completeness we have however to say that our analysis has been per- 
formed forcing the vector boson to decay into a lepton pair, as most of the experimental 
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studies do. For hadronic channels (i.e. W —>■ qq'), also the decay products are allowed 
to emit gluons and to give rise to a parton shower that is still described in the lead- 
ing soft /coUinear approximation by HERWIG. The implementation of matrix-element 
corrections to hadronic W decays is straightforward as they are very similar to the cor- 
rections to the process Z qq that are discussed in and is in progress. It is also 



worth remarking that the method applied in this paper to implement matrix-element 
corrections to the initial-state radiation in W and Z production can be extended to a 
wide range of processes that are relevant for the phenomenology of hadron colliders. 
Among these, the inclusion of matrix-element corrections to simulations of heavy quark 
production and particularly of top production in pp or pp interactions is expected to 
have a marked impact on the top mass reconstruction and many observables that are 
relevant for heavy quark phenomenology. This work is also in progress. 
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